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Introduction. 

§ 1. A consideration of the processes of cosmogony demands an extensive knowledge 
of the behaviour of rotating astronomical matter. What knowledge we have is 
based upon the researches of Maclaurin, Jacobi, Poincare, and Darwin. These 
researches refer solely to matter which is perfectly homogeneous and incompressible, 
although it is, of course, known that the primordial astronomical matter must be far 
from homogeneous and probably highly compressible as well. The question of how 
far we are justified in attributing to real matter the behaviour which is found to be 
true for incompressible and homogeneous matter is obviously one of great importance. 

§ 2. There are d priori reasons for expecting that there will be wide differences 
between the two cases. Consider first a sphere of homogeneous incompressible 
matter devoid of rotation. This will be stable if every small displacement increases 
(or, at least, does not decrease) its potential energy. The sphere has a number of 
independent possible small displacements which can be measured by the number of 
harmonics which can be represented on its surface. The spherical configuration is 
known to be stable because it can be shown that every one of these displacements 
increases the potential energy. 

Contrast this case with the corresponding one in which the matter is compressible. 
The number of possible small displacements in this latter case is measured by the 
sum of the numbers of harmonics which can be represented on all the spherical 
surfaces inside the sphere. Let R be the radius of the outer surface ; let r, r', r", . . . 
be the radii of all the spheres which can be drawn inside this outer sphere, 
and let r n , r f m r ff n ... R„ be the number of independent harmonics which can be 
represented on these spheres. To prove that the sphere is stable it is now necessary 
to prove that every one of the r n -\-r' n -\-r" n + ...R n possible displacements increases 
the potential energy. If we argue by analogy from the case of an incompressible 
sphere we are, in effect, merely considering R n of these displacements and neglecting 
the much greater number f n 4V n + r" n +. Furthermore, in these neglected dis- 
placements, the nature of the displacement is essentially different from that in the 
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R n displacements, so that there appears to be no justification at all for an argument 
from analogy. - 

In each of the neglected displacements, the change in the potential energy will 
consist of two terms. There will be a change in the elastic energy of the 
compressible material, and this can be easily shown to involve an increase in the 
potential energy. There will, in addition, be a change in the gravitational energy, 
and this can be shown # to involve invariably a decrease in the energy. If W, E, G 
denote the total, the elastic, and the gravitational potential energies, 

sw - se+sg, 

in which SG is invariably negative. The condition for stability is that for every 
possible displacement SE shall be numerically greater than §G. 

It might naturally be thought that by considering a system in which the matter 
was, so to speak, very gravitational or very little elastic we could have SE small or 
SG great, and so should have instability of the spherical configuration. But it must 
be remembered that the gravitation and the elasticity of the matter are not 
independently at our disposal. The action of the gravitational forces tends to 
concentrate the matter and so involves that the elasticity becomes large in the 
equilibrium configuration. If we consider a system in which the elasticity is 
artificially kept small, as, for instance, by adding an additional repulsive field of force 
to annul, or partially annul, the gravitational field, we can easily construct systems 
for which a spherical configuration is unstable, t but, short of this, it appears to be a 
general law that the elastic and gravitational agencies must march together in such a 
way that SE is always numerically greater than $G& so that every natural spherical 
system is stable. 

The nearest approach in nature to the artificial repulsive field imagined above is 
found in the influence of rotation. This influence may be represented by the super- 
position of the usual repulsive field of centrifugal force of potential —^w 2 (x 2 + y 2 ). 
The field is not spherical, and so the figures of equilibrium obtained under its influence 
cannot be spherical. But it can be regarded as made up of a spherical part —%w 2 r 2 
and a superposed harmonic disturbance |- w 2 ¥ 2 r 2 . The first term is certainly a 
spherical repulsive field, and will, of course, tend to annul the concentrating influence 
of gravitation. The problem which requires study is that of how far, or in what 
circumstances, the presence of rotation can disturb the otherwise general law that SE 
is always greater than §G. 

The problem is one of enormous complexity and great generality. It will hardly 
be expected that the present paper will contain anything approaching a general 

* Of. J. H. Jeans, "The Stability of a Spherical Nebula," 'Phil. Trans.,' A, vol. 199, p. 1. 
t Cf. J. H. Jeans, "On the Vibrations and Stability of a Gravitating Planet," 'Phil. Trans./ A, 
vol. 201, p. 157. 

J Cf. below, §§11, 22. 
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solution, and it may as well be stated at once that it does not. All I have been able, 
to do is to grope after general principles by solving a problem here and a problem 
there as seemed needful to illuminate a possible path towards a general theory, and 
the present paper is confined to a very few of the special problems I have considered, 
but I have selected those which seemed to have most bearing on the general question 
in hand. 

Medium in which the Pressure is a Function of the Density. 

§ 3. In the most general astronomical medium the pressure is, of course, not a 
function of the density. The relation between pressure and density varies from 
point to point, partly on account of inequalities of temperature and partly on account 
of variations of chemical constitution. But no general theory can be expected to 
apply to the most general heterogeneous mass of matter possible, and before any 
general theory can be deduced we must have material from which to deduce it. 

§ 4. The simple system from which we shall start will be a system in which the 
matter is homogeneous as regards its properties, so that at all points the pressure 
and density will be connected by the same relation. It will be seen later (§ 15) how 
it is possible, in at least one important respect, to escape from this limitation. 

For the present we assume the pressure and density to be connected by the 
relation 

at every point. We take the centre of gravity of the rotating mass to be the origin, 
and the axis of rotation to be the axis of z. The equations of equilibrium are 

ox ax 

dp 3 V ■ 2 

dy dy 

dp _ 3V 

dz " p a7 5 

in which V is the potential of the whole gravitational field of force, In virtue of 
relation (l), these equations have the common integral 



dp 



= V + ^ 2 (» 2 +^)+C (2) 



or 

<f>( p ) = Vh-^ 2 w 2 + C, . (3) 

in which w 2 stands for x 2 + y 2 y and <p(p) for I— ^, which is by hypothesis a function of 

j P 

3 N 2 
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p only. There is further the relation of Poisson, 

V 2 V = — 4tT/>, (4) 

so that on operating on (3) with V 2 we obtain 

V 2 (j) (p) + 4tt/) = 2w\ . ......... (5) 

the differential equation which must be satisfied by p in any configuration of 
equilibrium under a rotation w. 

§ 5. In general a solution of equation (5) will involve negative and zero values of p. 
In the physical problem p will be limited as to values, and this limitation will 
determine the physical boundary of the rotating mass. 

Let V m denote the gravitational potential at any point in space of the finite mass 

determined in this way. We have found a configuration of equilibrium under a 

potential V, the potential of the mass is Y m9 so that for equilibrium we require an 

additional field of potential V—V w . We can say that the configuration found will 

be a true configuration of equilibrium under an external field of force of potential V 

such that 

V = V 4-V (6) 

And, inasmuch as V 2 V m = —4^ = V 2 V, it is clear that V 2 V = 0, so that the 
external field has poles only at the origin or at infinity. The condition that any 
solution shall lead to a configuration of equilibrium for a mass rotating free from 
external influence is, of course, V = 0. 

§ 6. The simplest solution of equation (5) is obviously that in which p is a function 
of r only, but it is clear from (3) that this cannot give a free solution except when 
w = 0. 

§ 7. The next simplest form of solution is that in which p is a function of z and vs 
only, and this can give a free solution. It includes, of course, as a particular case 
the system of Maclaurin spheroids. For this class of solutions every section at right 
angles to the axis of z is circular, and in any such section the lines of equal density 
are circles. The density at any point is of the form p —f(m, z). 

Let denote colatitude measured from Oz, and let ^ be azimuth measured from 
the plane of xz. The most general configuration which can be obtained by displace- 
ment of that just considered will have a law of density of the form 



oo 



p=f {VS, z) + Z /, (tar, Z) COS S\lr. 

1 

It is easily seen that the separate cosine terms lead to independent displacements, 
and we shall for the moment only consider the displacement of the first order, for 
which the law of density is 

P =/o(<*> z)+/i(<*> ^)cos^, (7) 

where yi(isT, z) is a small quantity of the first order. 
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The boundary being a surface of constant pressure must also be a surface of 
constant density, say cr. The equation of the boundary is accordingly 

/ (m, z) +f x (sx, z) cos yf, = <7 (8) 

The whole mass inside this boundary may be regarded as composed of coaxial rings 
of matter as follows. Inside the figure of revolution f (m, z) = <r, we suppose there 
to be a series of rings of density given by (7), while the surface inequality can be 
regarded as represented by the presence of rings on this figure of revolution of 
density proportional to cos \fr. 

On integration the potential Y m at any external point is seen to be of the form 

where xo> Xi are functions of w and z only. 

Suppose now that the surface is so nearly spherical that spherical harmonic analysis 
may be used with reference to it, then, since Y m is a solution of Laplace's equation at 
all external points, and is also of the form (9), it must be of the form 

v„ = Ao + £ ( AA p; ( M ) cos ^ (io) 



T l \<y* s ~^~ 



where /ul = cos 0, and P/ (jul) is the usual tesseral harmonic — P s (jul). Moreover, since 

the centre of gravity of the mass is supposed to coincide with the origin, A 1 must 
vanish. 

We have, from equation (3), if V = Y mi 

= <j> {f (m, z)} +J[ (tjt, z) cos yfr <p f {/ (m 9 z)} -^V (l -P a ) -C 

at all internal points. Equating these two expressions, we must have at the 
boundary 



r 

CO 



/ 1 ( Wj ,).^{/ ( W) ,)} = s(A)p s i( M ) ! 



2 \r s+ - 



or, neglecting small quantities of the second order, 



/(w, *)*'(*) =Z(;^) P. 1 (m) (11) 



2 \r s+ 



Hence either J\ (m, z) vanishes at the boundary or is of at least the second order of 
harmonics. 

It follows that if there can be a configuration of equilibrium which differs from the 
configuration of revolution p =f (m, z\ by a displacement proportional to the first 
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harmonic, this configuration must be one in which fi(vt, z) vanishes at the boundary, 
so that the boundary must be a figure of revolution about the axis of rotation. 

It now follows from (3) that V must be a function of m and z only on the 
boundary, and hence also (since V is harmonic) at all external points. It follows 

that - — , and hence also ~- . are functions of tar, z only at the boundary. Whence 

again, by equations (4) and (5), it follows that ^-y and ~~ are functions only of ns 

and z at the boundary. And, by successive differentiation of equations (4) and (5), 
it is seen that all the differential coefficients of V and p are functions only of w and z 
at the boundary. 

It can be seen from this # that the configuration must be one of revolution through- 
out. In other words, there can be no configuration of equilibrium which differs from 
the configuration of revolution by first harmonic terms only. 

Laplace's Laiv. 

§8.1 have not found that any progress worth recording can be made with the 
general relation p = f(p), so that progress can only be hoped for by examining special 
oases. 

The case that suggests itself as most important is that of the gas law p ~ Kp, 
satisfied in a perfectly gaseous nebula at uniform temperature. The difficulty is that 
such a nebula extends to infinity in all directions, and so cannot rotate as a rigid 
body. Or rather, when it is caused to rotate, it throws off its equatorial portions 
and the remainder rotates in the shape of an elongated spindle of infinite length. In 
this connection I have worked out the purely two-dimensional problem of a rotating 
gaseous cylinder of infinite length. The results are too long to be worth printing ; 
it will, perhaps, suffice to record that the analysis bears out in full the conclusions 
arrived at in this paper. 

The law which is most amenable to mathematical treatment is Laplace's law 

or, as it is more convenient to write it, 

J> = ^ (p'-V), (12) 



K 2 



in which c, p, /c, and <r are constants, <r being the value of the density at the free 

* I have not succeeded in obtaining a rigorous proof of this. It might be objected that 
nothing in the above argument precludes first harmonic terms proportional to such a function as 

i 
e ~~' f(w,z) i where /(G7, z) — is the equation of the boundary. The pure mathematician may not, 

although the astronomer will, be influenced by the consideration that such functions never occur in natural 
problems. If such a function did occur, it would involve an extremely fantastic relation between p and p. 
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surface. This law has the merit that the case of an incompressible fluid is covered 
by the special value c = oo or k = 0, the density now having the value or throughout. 

There is the & priori objection to the law that its form precludes first harmonic 
displacements (cf. below, § 11). This objection would be fatal were it not that we 
have seen that first harmonic displacements are in any case of no importance. This 
being so, the objection falls to the ground, and I have thought it worth working out 
this law as far as possible . 

§ 9. Using the relation (12), we have in place of the more general equations (3), (4), 
and (5), the particular equations 

i^ p = V+>W*+C, (13) 

K 

V 2 V = -4tt/), (14) 



~ V 2 p + 47T/0 = 2w 2 . . (15) 



2 
K 



On putting 



this last equation reduces to 



w 2 



2 



IT 



~~ ~~ /v' * * * * * * * * * * * \ / 



(V 2 + k 2 ) x = 0. . (17) 

No Rotation. 
§ 10. When there is no rotation w = 0, x = P an( i the equation becomes 

The general solution is 

P = 2A w r~ 1/2 J w+ i /2 (/cr)S w , . (18) 

while the particular solution giving a spherical boundary is 

p ~ A r"" 1/2 Ji/ Ur) = 7== — , (19) 

the last being, of course, the well-known solution which occurs in Laplace's theory 
of the figure of the earth. It will now be shown that this configuration is stable for 
all displacements. 

§11. Let r = a be the free surface corresponding to the simple solution (19). 
Consider an adjacent solution 

P = A <) r- 1/ »Ji /s (/^) + A n r- ,/a J n+ i / ,(icr)S n , ...... (20) 

and let the corresponding free surface be 

r z=z ct ~t~ UK$ n * ........... v^^ / 
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On substituting this value for r in (20), neglecting squares of b and equating 
corresponding harmonic terms, we obtain 

<T = A C&~ Ji/ 2 \K(Xf) = — y= — , ....... (22) 

7 

A n a~ k J n+ i k (Ka) = — Aq6 -=- {a~ 1/2 Ji/ 2 (/ca)} = /cA 6a~ 1/2 Js /2 (ko), . . . (23) 
whence 

By integration, the value of Y m at a point on the sphere r = a is found to be 



a Jo ° ,2V ; (2w + l)aJ 







a" " T,2V ' 2w+l 



_ 47rA s; j / \ 4ttA„ M+ '/ 2 T ./^fl 



+ (2n+l)a B+1 " <tJ./,M 

while the value of V, as given by equation (13), is 

V = i^s «-'/. J 1/2 (,<«) + i^ a- 1 '' J, +v , M S re + cons. 

If we put 

V-V m = v = v +v n S„ 

we obtain, after some reduction, 

In general, this gives the value of A n which determines the tide raised by a field of 

-J S n proportional to S n . We notice that when n = 1, v n = or A n = oo 

independently of the values of /c, a. This merely expresses the obvious fact that 
there can be no equilibrium at all so long as the fluid is acted on by a force 
proportional to a harmonic of the first order. 

If it is possible for there to be a configuration of equilibrium when v n = 0, other 
than that given by A n = 0, this configuration will of course determine a point of 
bifurcation in the series of symmetrical configurations. The points of bifurcation are 
accordingly given by v n = -0, or by 

J w -i/ 3 (*fl) _ Jy» (*a) (oa\ 

-_. / \ ■— t ' / \ • * * • * * * * * * \ / 

J B + 1/2 (/CC*) J !/,(««) 

For brevity in printing, let us introduce the function u n defined by 

«» = H1 = -f:i Io ^W" ,+,,,j '-''>W}' • • • W 
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Near kg = 0, u n - 



2n + l 
the difference-equation 



; the value of u Y is 



K(X, 



cot *a, and successive ^'s satisfy 



Uft+i — 



2n + l 



Ka 



1 

u 



(28) 



n 



With the help of these properties, it is easy to draw approximate graphs of the 
curves y = u n . Such a graph, for values of ica up to the first zero of u x (jca = 4*49) is 
represented in fig. 1, in which the vertical scale is 2-g- times the horizontal scale. 




e27— O 



Fig. 1. Graphs of the functions u n . 



In terms of these functions, the points of bifurcation are given by u n = u v It is at 
once evident that there is no root of this equation for values of Ka less than 7r, and 
therefore (c/1 equation (19)) no point of bifurcation at all so long as p is restricted to 
being always positive. It follows that the spherical configuration is stable for all 
displacements. 

Small Rotation. 

§ 12. When the fluid experiences a slight rotation w 9 the spherical configuration is 
of course slightly flattened. The appropriate solution of equation (17) is 



X 



w 



i7T 



= A r- 1/2 Ji/. 2 (kt) + A 2 r- 1/2 J 5/2 (kt) P 8 , 



(29) 
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where P 2 is the second zonal harmonic about the axis of rotation as 6 = 0. Assuming 
the free surface to be 

the equations analogous to (22), (23), and (24) are found to be 

a-— — - = A a~ ll2 Jij 2 (Ka), . . (31) 

Zj7T 

A 2 a~ i! ' 2 J * i 2 (Ka) = /cA &a~ 1/2 J 3 / 2 (/ca), (32) 

&o- = A&-''' f^ ( ^] ^ M (33) 

\ Z7roy 

Let v be given by 

v = y+i w W-V w = V-V W +^V(1-P 2 ), 

then, instead of equation. (25), we have 

v = ilj^l ij 3k ( Ka )- J'fcM Jv-M \ p 2; (34) 

in which constant terms are omitted, and the value is taken on the sphere r = a. 
For a configuration of equilibrium under no external field of force we must have 
V = V m , and therefore v in equation (34) equal to — -Jw 2 a 2 P a . Neglecting squares 

of iv 2 , and therefore omitting the factor 1 — - — in the denominator, the equation 

becomes 

|^V = — A 2 a !i J ^(ko) (u 1 —u 2 ) > (35) 

Ok 

giving A 3 in terms of w 2 , when w 2 is small. It will be readily verified that this 
equation is identical with that obtained by Thomson and Tait ('Nat. Phil./ §824, 
equation (14)). 

§ 13. We next examine the solution 

x = / ,_^. = A r-^J 1/2 (/cr) + A 2 r- 1/ a 5/2 («r)P 2 + A B r- 1/ a B+v . 2 ( K r)S, i , . . (36) 

Air 

which is appropriate to a mass of fluid having a rotation w given by equation (35), 
and acted on by a field of force of potential v n $ n . By analysis exactly similar to 
that just given, we obtain at r = a 

__ 4:7rA n a ! * f T , v J w+ i/,(yq) Ji/,(*:a) "] / q7 x 

^""(2^+1) J J -vA^-7- ^x L OT 
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This gives A n for the general tide raised by the field v n $ n . The condition for a 
point of bifurcation is v n = 0, or 

^ = (1-—)^ (38) 

Thus the points of bifurcation, if any, are still determined by the intersections of 
the graphs in fig. 1, except that the graph of u x must be supposed decreased vertically 

in the ratio 1 — - — to 1. 

27TCT 

We may, if we please, imagine that we start with very small rotation, and allow the 
rotation progressively to increase, this increase being accompanied in imagination by a 
greater and greater flattening of the graph of u v 

It is clear that under all circumstances the curve which will first be intersected by 
the flattened graph of u x will be the graph of u 2 . It is further clear that the 
requisite value of iv 2 is least when ica = 0, and progressively increases as tea increases, 
at any rate up to tea = ir. 

This means that in the first place the circular vibration will invariably become 
unstable through a vibration proportional to a second harmonic, so that the first point 
of bifurcation reached will be one such that the spheroidal form gives place to an 
ellipsoidal form. If the rotation is so small that the problem may be treated as a 
statical one, there will be no question as to there being an actual exchange of 
stabilities at the point of bifurcation, for clearly v n changes sign at this point. Thus 
for rotation greater than that at the point of bifurcation, the spheroidal form will be 
definitely unstable, and the ellipsoidal form definitely stable, at least until the next 
point of bifurcation is reached. 

Our result shows, in the second place, that the masses which become ellipsoidal for 
the smallest values of w 2 are those for which ica is smallest. To put it briefly, the 
mass which is most unstable when it begins to rotate is the incompressible mass — a 
somewhat unexpected result. 

For any value of ica, the value which w 2 must have for the spheroidal form to 
become unstable is \cf. equation (38)) 



w 2 %u 



i-^ = f> (39) 

and when /ca = 0, the value of uju x = § (cf. § 4). 

Thus our equations would make the spheroidal mass of incompressible fluid first 

become unstable when = "400, but these equations have only been obtained on 

27TCT 
VI) 2 • • • 

the supposition that — — is so small that its squares may be neglected, a supposition 

which is now seen d posteriori to be hardly admissible. Probably the results 
obtained are qualitatively true, but quantitatively unreliable. In point of fact the 

3 o 2 
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first point of bifurcation for an incompressible mass, instead of being given by 



= '400, is known to be given by the widely different value = '1871. 



Our analysis has nevertheless proved rigorously the point which is really most 
important, namely, that there can be no point of bifurcation at all for quite small 

values of - — . At the same time, since the question of when and how a rotating mass 

27TCT 

first becomes unstable is one of considerable importance, I have attempted to obtain 
a more reliable investigation than the preceding. I have found that the accuracy is 
not greatly improved by taking the analysis as far as squares of w 2 l2ir(T, while the 
labour of working with a general power series would be appalling. I have, therefore, 
reluctantly been compelled to give up hopes of carrying the rigorous solution of the 
problem further in this direction, but have thought it worth while to examine the 
analogous problem for rotating cylindrical masses. All the essential physical features 
of the natural three-dimensional problem appear to be reproduced in the simpler 
cylindrical problem, so that it seems legitimate to hope that an argument by analogy 
may not lead to entirely erroneous result. 

Cylindrical Masses in Rotation. 

§14. The fundamental equations are, of course, the first two of the equations 
already written down in § 3. The third equation does not occur, since — = 0. The 

equations have, as before, the integral (13) leading to the differential equation (15) 
for p. 

The most general solution possible will be 

p — y 2 A n J n (Kr) cos (nO—e), (40) 

27T 

in which r now stands for <sJ(p(?-\-y 2 \ No matter how great the rotation, there is 
always a special circular solution 

p = — +A J (kv), (41) 

2tt 

this being analogous to the spheroidal figures of equilibrium investigated in § 12. 
Let us examine the deformed solution 

p = jr- + A J (kt) + A W J W (kt) cos n6, (42) 

/ w 2 

in which A n is supposed small, but there are no restrictions on the value of — . If 
the free surface p = a- is supposed given by (cf. equations (21) and (30)) 

r = a + b cos n6, 
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then, as in equations (31) and (32), a and b must satisfy 



w 



or = -- + A J (ko). 

Air 



(43) 



whence 



A n J n (/ca) = — A 6/c J ; (/ca) — A bicJ 1 (ko) 



* 1- 



w 



i7T, 



J x (*a) 



(44) 



The potential of the mass, V m , can be regarded as arising from a distribution of 
density p inside the cylinder r = a, together with a surface density bar cos nO spread 
over the surface of the cylinder. 

The first part of the potential, evaluated at R, 6, is 



C- [log {r 2 + R 2 -2rR cos (6-0)}] 



v? 



i7T 



+ A J (kt) + A w J w (kt) cos n6 



r dr d6 



= C-2 



00 



/y» a 



sBf 



w 



log R— 2 -=5- cos 5 (6—Q) — + A J (*r) + A n J w (*r) cos n# 



(7T 



r dr d0 



— A w 



'a 



2tt r^ 
w R 



J w (kt) r dr cos ?i0 + terms independent of 



2 



7T 



*n 



A-^Jw+i ( lca ) ~f^r cos ^0+ terms independent of 0. 



The potential of the surface distribution is 



so that, at r = a, 



27rcrba n+: 
nR n 



cos nQ, 



V m = i A n JL +1 (*a) + — - [ cos n0+ terms independent of 0, 

{ m n ) 



while, by equation (5), 

V = -~ \ 7T- + Ao J (/cr) + A n J w (/cr) cos n# f — \w 2 r\ 

K [ZtT J 

If, as before, we express the tide-generating potential V— Y m in the form 
v + v n cos n6, we obtain for the value of v n9 at R = a, 



v 



w 



a 1 4-7T T / \ 2-7ra T / \1 27rbara 

k ten j n 



kU 



-n 



2-7ra a f T / \ Jo (k<%) J w (ko) 

\ Air (J J 



< 



* (45) 



It will be seen that this equation is exactly analogous to the former equation (37), 
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but with the important difference that the present equation is true for all values of 
w 2 , without limit. The points of bifurcation are given by v n = 0, or 



w 



u n _ij 2 (ka) = 1 - - — ) Ulrica), (46) 



iir<y I 



which again is exactly analogous to the former equation (38). The graphs of the 
functions % 2? % 2 , . . . will be found to lie as in fig. 2, and we may again imagine that 




ac=o 



z 3 

Fig, 2. Graphs of the functions ^ w+ i/ Q 



4 



points of bifurcation are sought by flattening the curve U\^ until it intersects the 
other curves. 

It is clear that, under all circumstances, the first curve to be intersected will be 
the curve % 2 , corresponding to a displacement proportional to cos 29. Thus, as 
before, when the circular form becomes unstable, it gives place to a form of elliptic 
cross-section, which is stable. Moreover, the smaller ko, is the lower the value of 
w 2 /27rcr for which the circular form becomes unstable. 

These results are true without any regard to the value of w\ so that they confirm 
the results stated, but not rigorously proved, in § 5. The numerical calculations 
which follow will make the matter clearer. 

If p denotes the mean density of the rotating mass, the total mass per unit length 
is given by 



pira 2 



»7T 



'a 

prdr = 2ir 



w 2 a 2 A a n J 1 (ko)\ 
- — + a r , 







:7T 



K 
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giving, on substitution from equation (43), 



w 



p = h 



Ji ( Ka ) 



W 



2lT KG J (kO) \ 27T/ 2-7T 



w 2 . r 1 ( j 2 (ko) \ ( w 2 



= HL + ( 



\ 



J (ica)/ \ 2 



7Ty 



whence, for the ratio of p to <r, we have the general formula 



£ = 1+ J 2 (ica) U ™* \ t 
a- J (ica)\ 2xflr/" 

For the particular configuration which occurs at a point of bifurcation, 



1- 



W 2 __ J (ko) J 2 (/ca) 



i7T0" 



so that 



J x 2 (kO) 



J 2 (^) 



whence we obtain 



£ - 1 + 

a- (Jj (/cCfc). 



i + A 2 , 



w 



2 



^k-U*h 



2irp -Uii 3 (l+u\y 



In the following table I have calculated the values of w 2 /2ttct and of w 2 /2wp for 
which cylinders of different radii (a) and compressibility (k) first become elliptical in 
cross section : — 



Ka. 


V 


3 /2 


w 2 

27rcr' 


w 2 
2irp 



i 

3 

¥ 
1 

2 
2 • 4048 

3 
3-8317 


o • 0000 
0-1261 
0-2582 
0-4040 
0-5751 
2-575 

00 

-1-304 
0-000 


0-000 
0-063 
0-126 
0-192 
0-261 
0-612 
0-829 
1-433 

CO 


0-500 
0-503 
0-510 
0-525 
0*546 
0-762 
1-000 
2-099 

CO 


0-500 
0-500 
0-502 
0-506 
0-511 
0-554 
0-593 
0-687 
1-000 



It will be seen that the general result is fully confirmed, that incompressible 
masses are the first to become unstable, and that the more compressible the mass is, 
the greater is the rotation required for it to depart from a symmetrical configuration. 



Rotating nearly Spherical Mass with High Internal Temperature. 

§ 15. We now leave the artificial two-dimensional problem and return to the real 
problem in three dimensions discussed in § 13. 

The coefficient tc was there assumed to have the same value throughout the mass, 
as of course it would if the matter were homogeneous and of uniform temperature 
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throughout. But to represent natural astronomical conditions there is no question 
that ac ought to increase on passing from the centre to the surface, thus representing 
a mass in which the temperature is highest inside and falls towards the surface. 

We are in this way led to study the question of stability when k is a function of r. 
It would be difficult to say precisely what function ought to be chosen if we were 
trying to represent natural conditions as faithfully as possible. It appears, however, 
that no continuous function will lead to equations which admit of integration. The 
only case which appears to be soluble is that in which the matter, before rotation, 
may be treated as if formed of a series of different layers, each being homogeneous 
and at a uniform temperature in itself, but the temperature varying from layer to 
layer. To represent this we take different values of k in the different layers, k being 
smallest in the interior. 

There is no limit to the number of layers which can be treated analytically, but 
the assumption of a great number of layers naturally leads to highly complicated 
formulae which are capable of conveying their meaning only after laborious numerical 
calculations. Both in order to obtain comprehensible results and to simplify the 
argument, the layers will, in what follows, be supposed to be only two in number. 
They may conveniently be referred to as the core and the crust. It will be found 
possible to generalize the results obtained so as to apply to any number of layers. 

§ 16. We accordingly suppose that there is an interior core of radius a, in which 
the coefficient of compressibility has the uniform value /c, and that outside this is the 
crust of external radius c, in which the coefficient is k . It is again necessary to 
suppose the rotation to be so small that w 2 may be neglected. 

As in § 3, the density p must satisfy 

(V 2 + k 2 )L~^)^0 (47) 

throughout the core, and the same equation with the appropriate value of k throughout 
the crust. The most general solution of equation (47) is 

P = %- + 2 {A»r- 1/ 'J n+ . / ,(icr)+B ) ,r- v «J_ 6l+ i W (Kr)}S I , (48) 

In the former problem all the terms in B n could be omitted because p had to be 
finite at the origin. In discussing the solution for the crust these terms must be 
retained. The solution can, however, be put In a more concise form. 

Let the constants A*, B n be replaced by new constants G n9 6 n given by 

A n = C n cos n , B„ = C n sin Q M 
and let us introduce a function J n+ i/ 2 (x, 6) defined by 

J»+V. (^ e ) = J »+V* ( x ) cos e + J-0.+ 1 /*) ( x ) sin 0. ..... . (49) 
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Then the solution (48) may be replaced by 



w 



J-J 



7T 



4- 2 C tt r- ,/9 J tt+1/2 (kv, 6 n ) S n , ....... (50) 



which is formally analogous to (7). 

The following properties of the function J n+ i /3 (#, 0) may readily be verified, and 
will be required later : — - 

£{a^W . . . . . (51) 

d {ar ( » + WJ, +Vt (a, 0)} = -ark+WJ^fo -0), ..... (52) 



c£,x 



2n+l 



There is a ready rule for writing down the values of these functions. In the first 
place, we have 

Tf \ Dill Ju T / \ COS «Xy 



so that 

T , . sin(aj + 0) 

J ,,„(#, 0) = 7T= ^ e 

Now let 96 (# + 0) be used to denote a general function made up of circular functions 
of x + 6 and of algebraic functions of x. Then Ji/ S (#, 6) is of the form </>(% + 6), and 
any number of differentiations with respect to x, or of multiplications by powers of x 9 
will still leave it in this form. It follows from (52) that J 3/ , (x, — 0), J« /a (aj, 0), &c, 
will be of this form. Hence we have the general law 

J w+ i/ 3 (^, 0) = $ {^-i-(~~l) n 0} ? . ...-..,. (54) 

in which the functional form of <f> is at once given by 

(a?) = J* + i/ 2 0*0- 
For instance, 

J 5/ (#) = ( — — 1 ) sin as cos x> 

,2X } \x 2 J x 

so that 

•M*. 0) = ,(* + *) = ft-l)sin(, + 0)- 3 -cos(x + e).' 

\*JU / X 

§ 17. We proceed to carry out analysis similar to that of § 13. Suppose that under 
a tide-generating potential v n 8 n} and a rotation w 9 the core assumes a configuration 
such that its boundary r = a becomes deformed into 

/ == l$ ~f- OhJ n ~T" fj Jl g ? . . , . . » , , , , ( «J ) 

VOL, OCXIII, — A, 3 P 
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while the outer surface becomes 

r = e + cE3 n H dP 2 , ........... (56) 

Let us suppose the densities in the two layers to be 

2 

/. = ^-+A r- 1/ 'J Vt (ic?-)+A,r- ,/ 'J 11+v ,(icJ-)S l , + A^r- ,/ 'J., a ( K »-)P 3 • • • (57) 
in the core, and 

/» = ^ Cor-'^J,/, («r, a) +C„7-V.J B+1/a (*,-, ,8) S.+ C^-'M.,, (, t -r, y ) P 3 . , (58) 

in the crust. The boundary between the two layers must clearly be an equipotential, 
and therefore a surface of constant pressure and density. Let o~, a be the densities 
at this boundary in the core and crust respectively. 

On replacing r by a + bS n + /3~P 2 in equations (57) and (58), the values of P must 
become a- and <r respectively. This leads to the relations 



w 2 



G" — 



2 



7T 



J±.qC(/ 2 uif, 2 \KCif) 9 , , . . . . , . . . (59) 



iv 2 



,/_ _ — C a 1/-2 Ji; 9 (//a 5 a), .......... (60) 

\ 27T/ /cJ3/ 2 (/ca) 

/_ ^!A = C a _1/ - 2 ^~ n+l/ * ^ a " ^ ^ ^ a " a ^ (62) 

From, similar analysis applied to the outer boundary, if <x is now the density at 
this boundary, 

— C C - ' 2 Ji/ 3 (/c'c, a), .......... (63) 



""o 



2 



dL-^) = c n c- i i*± * +i i4i c fi Jii 'H c ' J ^ (64) 

\ 27r/ k Js/ 2 (jcC, --a) 

Similar equations, of course, connect the coefficients which depend on the rotation. 
The value of Y m at a point on the sphere r = c can now be written down, as in 
§11, and is found to be 



4?r 



V = 



— — + -~ {c 3/a Js/^/c'c, — a) — a 3/ %L/., (//a, —a)} + — a< v \L/,(/ca) 

ZlY OK s " " /C 

CI 



h (2rn-l)«- +1 "' "'• " l ; 
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in which the rotational terms proportional to the second harmonic are omitted and 
U„ is given by 

K K 

+ da- C n+2 + b(a—a')a n+2 (66) 

The value of V at r = c is, from equation (13), 

v = -%\z-+ G oc~ 1 '* J>/. ( KG > a ) + G « c " 1/a J »+v. (^ Z 3 ) s » r -i^ 2 (* 2 +2/ 2 ) + °o ns -> • ( 6 7) 

/c LZtt J 

whence, evaluating V— Y m and picking out the coefficient of S n , we find as the value 
of v n at r = e, 

*>» = ^5 n C- ,/a J n+ i /s (/cC, /3) - 72^ + ^\ c n + i U ** • • • • • ( 68 ) 

As before, the points of bifurcation, if any, are given by v n == 0. 

§ 18. It is now necessary to consider the boundary conditions which must be 
satisfied at the junction of the two layers. The condition of continuity of material, 
i.e., that the inner surface of the crust shall coincide with the outer surface of the 
core, has been expressed in equation (55), b and ft being the same for both core and 
crust. There is an equation of continuity of pressure expressed by 



2 2/22 

Q" ~-Q"oo __ Q" ~-Q"oo 

' "~" /2 

K K 



(69) 



which o- 00 is now used to represent the density associated with zero pressure. 
Finally, there is a condition of continuity of normal force and this requires careful 
discussion. 

Let M x denote the mass of matter actually forming the core and let Y 1 denote its 
potential at any point outside the core. Let M 2 denote the mass which would 
replace the core if the solution (58) for the crust were extended to the centre and let 
V 2 denote the potential of this mass at any external point. It will be noticed that if 
solution (58) were extended to the origin, it would give an infinite density p at the 
origin and also an infinite value of V in virtue of equation (5). On the other hand, 
it is readily found, by direct integration, that V 2 the potential of the mass M 2 , is 
finite at every point, including the origin. It follows that V can only be the 
potential of this imaginary arrangement of matter when it is acted on by certain 
external forces of which the potential becomes infinite at the origin. Let V 3 represent 
the potential of these forces. The value of V 3 is readily found, for it must satisfy 

V 2 V 3 = 0, and must coincide with V or with — w -~ to within an additive constant at 

K 

the origin. Thus V 3 is the limit of the right-hand side of equation (67) when r = 

3 P 2 
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r replacing c. This is found to be 



V, 



/c 



:7T 
12 



" C (Jtt)^ sin a G n (jir) lk l. 3.5... (2n-l) (~l) n sin )3 g 



n 



The condition now to be satisfied is clearly that 



at all points on the boundary r = a + &S n . This requires that Y 1 — V 2 — V 3 shall vanish, 
to within a constant, at all points outside this boundary, and therefore, in particular, 
at r — c. It will be readily seen that 

v 1 v 2 — v m \ v mja = 

while V 3 is exactly the value of the terms in Y m that involve a, when a is put equal 
to zero. The conditions sought are, therefore, simply that all the terms in a which 
occur in V m shall vanish at every point of the sphere. 

§ 19. We may now equate the coefficients of the separate harmonics, and obtain 



~Y J 3 / 2 (//a, — a) = — - Js/. 2 (ko) , . 



t • * * * 



• (70) 



%J n+ ,,,( K 'a,-/3) + &(,/-^W / ' 



/C 



27T/ 



A 



/c 



J n+ »/ 3 (*«) + & K 






• • 



(71) 



On account of the simplifications made possible by these equations, equation (68) 
may be put in the form 



'n 



n /2 

K 



4?T p 1/ r / „\ 4?T \C n +8 / y / / ~\ 7 »+2 I 



L \ ZwcrJ 



* * 



• ■• ■• 



o> 



The elimination of A and C from (59), (60), and (70) gives 



z w 2 \ Js/ 2 (/c 7 ^, —a) 
27T/ kSiiJjcO,) a) 



27T/ /C J i/ 2 (k6&) 



while similarly the elimination of A n and C n from (61), (62), and (71) gives 



w 



\ \ l J n+ > u (ic'a, -ft) J./.Qc'a, -*) \ = 



2tt/ ) J n+1/3 (*'a, /3) J Va (#r'a, a) J 



w 



2\ 



<T 



2 



1 + 



7T 



J w+8/g (ya) J» /a (ica) \ 
J w+ i/ 2 (/ca) J V2 (*a 



(72) 



,-^iH, ...... (73) 



. . (74) 



For brevity we introduce a function t& w (x, 0), a generalisation of the u n of § 11, the 
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new function being defined by 



u n (x, e) = -A- io g {x-*+*j ( x> 6 )} = J^fe -0) , . . . ■ ( 74a) 



Equations (73) and (74) now become 



/ w 2 w 2 

cr <r — 



7 % (/c'a, a) = U x (ko), (75) 



K K 

/ ..,2\ / „..2 



<r'- ~) {l+'ih (*'a, a) u n+1 (/c'a, /3)} = (a-- ~J {1+^ (*a) w n+1 (*a)}, . (76) 
while (72) becomes 

*, _ 47TC 1/g C n r , , ^ f 1 1 U h (k'C 9 -ft) "] / 77 N 



irar 



so that points of bifurcation are given by 

,2 \ 



w 



u n ( K 'c 9 -p) = 1-^-^(^,4 ....... (78) 



27ra" 



Again, if the rotation may be treated as small, there will invariably be a change of 
stability at these points of bifurcation, since vJG n changes sign on passing through 
one of them. 

§ 20. It is at once clear that the method can be extended to a mass consisting of 
any number of layers- — the only difficulty occurs in the numerical computations at 
the end. At each boundary between two consecutive layers there will be equations 
of continuity precisely similar to (75) and (76), while the final value of v n will be 
given by an equation exactly similar to (77), which it will be seen involves only 
quantities associated with the outer boundary. 

The procedure in any particular case will be to start, so to speak, with the 

innermost core of the system. Equation (75) is linear in cos a and sin a, so that 

tan a is uniquely determined. Leaving out of account systems in which the density 

is, in any part, negative, this will be found to be adequate to determine a uniquely. 

Equation (76) now becomes a linear equation in cos /3 and sin /3, from which /3 can be 

determined uniquely. In this way, passing from layer to layer, we can determine 

the various values of a, /5 for the different layers. Finally, the a's and /3's being 

known, equation (78) can be regarded either as an equation for w 2 or as an equation 

for c, i.e.) it can be regarded either as determining the highest rotation for which 

the symmetrical configuration is stable for a given value of c, or as determining the 

largest value of c for which the mass is stable under a given rotation. If the value of 

')jf} 2 . • ... 

— obtained by the first method is not small, the result will be inaccurate; if the 



27TCT, 







value for e obtained by the second method is so great that the density is in places 
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negative, the result will be of no interest except as proving stability for smaller 
values of c. 

§ 21. It may be well to take a general survey of the equations before giving special 
calculations. For simplicity we again consider two layers only, core and crust. From 
(75) and (76) it is clear that, when a = or k = k (involving o- — </), the values of a 
and j8 vanish. Broadly speaking, the more distinct the core is from the crust, the 
larger a and (3 are. Equation (78), of course, differs only from, the corresponding 
equation previously found, by the presence of the terms a and /3. The effect of these 
terms is seen on noticing that, in the notation already used, %i Y (kg, a) is of the form 
<j> {kg — a). Thus, to allow for the effect of the core on the term Ui(kg, a), we have to 
leave the algebraic part of the function unaltered, but to change all the trigono- 
metrical arguments from kg to kg — a. Speaking very broadly, the general effect on 
the graph of u x (cf fig. l) is a compromise between leaving the graph unaltered and 
moving it bodily a distance a along the axis. Similar statements apply to the graph 
of u n . Thus, while rotation as before is represented by flattening the graph of u x in 
fig. 1, the presence of a core is represented by a distortion of the graphs which may, 
with some truth, be thought of as bodily movements parallel to the axis. These 
bodily movements may cause new intersections between the graph u x and the other 
graphs, and the points of intersection will represent points of bifurcation at which 
the symmetrical configuration will become unstable. 

No Rotation. 

§22. The case that may properly be inspected first is that of no rotation. The 
equations reduce to 

— u x (//a, a) — - u x (ko) > (79) 

K K 

<t {1 + u x (//a, a) u n+1 (*'a, .0)} = a- {1 + Wi (tea) u n+l (ko)} 9 . . . (80) 

and, the equation for points of bifurcation, 

^^(/c / 6', — /3) = '^(kg, a) . (81) 

When n — 1, it is seen that /3 = —a is a solution of (80), and must therefore (§20) 
be the only solution. To verify that /3 = —a is a solution, replace /3 by — ain(80) 
and it becomes 

M i+ Jv 2 (^ a )j™"r + J V ,MJ' (82) 

which is seen to be identical with (79) (cf. equations (53) and (74a)). Equation (81) 
now reduces to an identity, so that every configuration is formally a point of 
bifurcation. The interpretation is, of course, the same as that of § 11, the displace- 
ment for which n = 1 is a rigid body displacement, and so requires no force to 
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maintain it. There is, of course, nothing of the nature of a change of stability, for 
vJG m instead of changing sign, remains permanently zero. The consideration of 
n = 1 is of no value except that it provides a check on the result of a rather involved 
series of analytical processes; 

When there is homogeneity between core and crust the non-rotating system has 
been found to be stable for all displacements. To examine whether this is altered by 
the presence of the crust, it is natural to test first the extreme case in which the 
difference between the core and crust is as great as possible. Let us make the core 
so hot that its density is zero, so that k has to be zero in order that the internal 
pressure may be maintained (of, equation (12)). 

Putting a- = ? equations (79) and (80) reduce to 



Uj (ic'a, a) = 0, u n+1 (//a, ft) — °°, . 



» * 



• * 



• (83) 



or, by equation (74a). 



J«/ 2 (//a, -a) = 0, J w+ i /2 (//a, £) = (), . 



8 # 



• (84) 



whence (equation (49)) a, ft are given by 

J _»/ a (k a) 



tan ft 



J„( n4 J/,)(/cV) 



» t » » 



(85) 



The values of a, ft corresponding to a few values of ko, are given below- 



KG. 











/ 










1 


- 12 


18 


2 


51 


9 


o 


100 


20 


4 


153 


13 


5 


- 207 


48 


6 


263 


12 


7 


319 


8 


9 


- 432 





12 


-- 602 


21 




fi- 



ll 



3. 





59 

15 7 

48 12 

91 54 

140 45 

192 21 

245 33 

354 51 

521 54 





2 

2 21 

16 39 

46 22 
86 9 



The case which is most favourable to the occurrence of points of bifurcation with 
positive values of p is when o- falls to zero at the outer boundary. Let us accordingly 
examine this case. We have (equation (63)) 



so that 



Ji/., (kc 9 a) = 0, 



K C 



7T — a. . 



. (86) 
. (87) 
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And in virtue of (86), the equation giving points of bifurcation (going back to 
equation (72)) is 

Vn . = 2^1 "7 J *-'I*\« C > ~~^> V 88 ) 

so that points of bifurcation of order n are given by 

J M _i/ 2 \7f — At, — p) = 0. , , . , , , . , (89) 

When n = 2, this becomes 

tan (tt-~ a + 0) — ir—a ; 
when n — 3, it is 

tan (tt— a— j8) = - — V — \V 

3— (tt— af 

On treating these equations numerically it is found that they can never be satisfied 
We conclude that the non-rotating mass is stable for all displacements, subject, of 
course, to the condition that the density shall be everywhere positive. 



Slow Rotation. 
§23. We consider next the stability of a rotating mass of the type under 



9, 



consideration, in which we are limited to — being small compared with the density of 

27T 

the main mass. If we suppose that o-, the density of the core at its outer boundary, is 
equal to — , we shall have a case— somewhat artificial of course — -in which the density 

of the core is very small compared with that of the crust, and in which the equations 
are not too complex to admit of treatment. 

We accordingly assume that a- = — , and the equations (75), (76), and (77) (or (72)) 

Zir 

reduce to the same equations as in the case of no rotation (equations (83)), Thus a, /3 
have the same values as before, being given by the table on p. 479. 

If we suppose that at the outer boundary of the crust the density falls to the small 

value <t = — , then the value of c, the radius of the outer boundary, is, as before, given 

by equations (86) or (87), and the value of v n is still given by equation (88). Thus the 
analysis is exactly the same as in the case of no rotation, and there are no points of 
bifurcation. 

It follows that, when a- does not have this special value assigned to it, the only 
hope of finding points of bifurcation rests upon the gravitational tendency to instability 
which arises from the presence of the small layer of crust in which p has a value less 

than — . Let us pass at once to the examination of the extreme case in which <r = 0, 

27T 
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Denoting as before the density of the crust at its inner surface (r = a) by <j\ and 
putting <r , the density at the outer surface (r = c) equal to zero, we have 



/ VJt 



C a 1/a Ji/ 2 (*'a, a) = cr 7 — -i 



2 



7T 



C c 1/s Ji/ 2 (kC 9 a) 



w 2 



2 



7T 



whence, on elimination of C, 



0? 



2lT<J 






c sin (/a + a) 



W 



a sin (/c + a) * 



. . . (90) 



Equation (72) still gives 



J„_l/ 2 (*'c, — /?) = 0, 



as the condition for points of bifurcation, and when n — 2 (the only case which 
appears to be worth examining), this reduces to 

tan (kc + $) = kg 9 .......... (91) 

in which ft is given from the table on p. 479. The procedure is to find kc from equation 
(91), and hence calculate the value of w 2 J2W from equation (90). The results for a 
few values of ko, are given in the table following (the last column is explained 
later) : — 



Kd, 


KC. 


cja. 


tV 2 j27rcr'. 


w*/2ir0. 





4-489 - 257 27 


00 


CO 


0-40 


1 


4-475 = 256 54 


4 


•475 


0-2221 


0-45 


2 


4-713 = 269 21 


2 


•356 


0-2269 


0-50 


3 


5-380 = 308 16 


1 


•793 


0-2159 


0-52 


4 


6-153 - 352 37 


1 


•538 


0-1821 


0-51 


5 


7*026 = 402 40 


1- 


405 


0-1570 


0-50 


6 


7-944 = 455 7 


1 


•324 


0-1365 


0*49 


7 


8-886 - 509 9 


1-269 


0-1198 


0-48 



The obvious remark must at once be made that probably all the values for w 2 J2ircr f 
are too large for results obtained by the neglect of w* to be accurate. But apart from 
absolute accuracy there is an obvious tendency for the value of w 2 J2W to fall off as kg 
increases — for lower values of cja the symmetrical configuration becomes unstable for 
lower and lower values of w 2 J2W. For kcc = 100, the value of w 2 /2W is 0*0104. 

§24. Against this, it must be noticed that the value of w 2 J2ira- f is of very slight 
importance ; what we are concerned with is the ratio of w 2 /2ir to the mean density of 
the whole mass. For a very rough calculation, we may assume the mean density of 
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the crust to be |-or ; , whence it follows that the mean density of the whole mass will be 
roughly equal to a density 6 defined by 

>3 „S 






e 

c 



and the value of w 2 /27r9 will be approximately the same as the quantity tv 2 /2irp which is 
computed from observations of binary stars. Values of %v 2 J2ir6 are given in the last 
column of the table on the preceding page ; the value 0*40 corresponding to kg = 
(no core) being inserted from the result of the previous analysis (§13). As before, 
the numbers are not numerically accurate, but their general trend may be expected to 
reveal the general trend of the true series of numbers. It at once appears that the 
values of w 2 f2irQ are surprisingly steady : there is certainly no rapid decrease in their 
amount as tea increases. 

Summary and Conclusion. 

§ 25. The problem we have had under consideration has been that of testing whether 
the behaviour of a rotating mass of compressible heterogeneous matter differs very 
widely from that of the incompressible homogeneous mass which has been studied by 
Maclaubin, Jacobi, Poincake, and Darwin. The result obtained can be summed 
up very briefly by saying that the ideal incompressible mass has been found to supply 
a surprisingly good model by which to study the behaviour of the more complicated 
systems found in astronomy. The problem especially under consideration has been 
that of determining the amount of rotation at which configurations of revolution 
(e.g.y spheroids) first become unstable. In so far as we have been able to examine the 
question, it appears probable that the compressible mass will behave, at least up to this 
point, in a manner almost exactly similar to the simpler incompressible mass, and 
results obtained for the latter will be nearly true, both qualitatively and quantitatively, 
for the former. The compressible mass, set into rotation, will apparently pass through 
a series of flattened configurations very similar to the Maclaurin spheroids ; it will 
then, for just about the same amount of rotation (as measured by w 2 fp), leave the 
symmetrical form and assume a form similar to the Jacobiah ellipsoids. Beyond this 
stage our analysis has not been able to deal with the problem. Indeed, strictly 
speaking, our analysis has hardly . been able to carry this far. A question of 
importance has been whether the quasi-spheroidal form for a compressible mass does 
not become unstable for a much smaller value of w 2 than the incompressible mass, and 
whether the instability does not set in in a different way. These questions we have 
been able to answer, with, I think, a very high degree of probability, in the negative. 
The whole matter is of necessity one of probability only, and not of certainty, for the 
general heterogeneous compressible mass is not amenable to analysis until a great 
number of simplifying assumptions have been made. 

It was first pointed out that a compressible mass has an infinitely greater number 
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of vibrations than an incompressible one, and as the mass is only stable when every 
vibration individually is stable, it might be thought that a compressible mass had more 
chance of being unstable — or would become unstable sooner — than the corresponding 
compressible mass. This has on the whole been found not to be the case, and on 
looking through the analysis the reason can be seen. 

A vibration in a compressible mass may be regarded loosely as a system of. waves ; 
the distance from one point of zero displacement to the next may be regarded as a sort 
of wave-length of the vibration. The stability or instability of a vibration depends on 
which is the greater — the gain in elastic energy or the loss in gravitational energy 
when the vibration takes place. But as between a vibration of great- wave-length and 
one of short wave-length there is this important distinction : for equal maximum 
amplitudes the gravitational disturbance caused by the disturbance of great wave- 
length is much greater than that caused by the disturbance of short wave-length, since 
the elements of the latter very largely neutralise one another. Thus the change in 
gravitational energy is enormously the greatest for disturbances of great wave-length, 
while it is easily seen that the changes in elastic energy are approximately the same. 
It follows that if the mass becomes unstable it will be through a vibration of the 
greatest possible wave-length, i.e., a wave-length about equal to the diameter of the 
mass. This general prediction is amply verified in the detailed problems that have 
been discussed. When we reflect that the vibrations of greatest wave-length are 
exactly those which ai*e common both to compressible and incompressible masses, we 
see readily that, in this respect at least, compressibility is likely to make but little 
difference. 

The vibrations of greatest wave-length are put in evidence, both in the compressible 
and incompressible mass, by the displacement of the surface. A vibration in which 
the displacement is proportional to a zonal harmonic P„ may be thought of as having 
a wave-length approximately equal to wa/n. In accordance with the principle that 
vibrations of great wave-lengths are most effective towards instability, we should 
expect the lowest values of n to give the vibrations which first become unstable, and 
this is, in fact, found to be the case. But here a very real distinction enters between 
the compressible and the incompressible mass. In the incompressible mass vibrations 
of order n = 1 are non-existent, the displacement being purely a rigid body displace- 
ment ; in the compressible mass vibrations of order n — 1 can certainly occur, and so 
might reasonably be expected to be the first to become unstable. 

It is in point of fact known that the incompressible mass becomes unstable through 
vibrations of orders 2, 3, ... in turn; it is found that the compressible mass also 
becomes unstable through vibrations of orders 2, 3, ... in turn, the vibrations of 
order 1 failing completely to produce instability. The reason for this apparent 
anomaly can, I think, be traced in the following way. In a displacement of order 1 
any spherical layer of particles will after displacement be spread uniformly over 
another sphere excentric to the first. The gravitational force produced by this 
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sphere of particles both before and after displacement is exactly nil at a point inside 
the sphere. Thus the gravitational field set up by a displacement of order 1 
neutralises itself in a way not contemplated in the general argument outlined above. 
Also the vibrations of order 1 and of greatest wave-length in the interior are not 
available, for they represent solely a rigid body displacement. 

The question of vibrations of order 1 is treated in §§ 3-7 ; it is shown that they 
may be disregarded, and we pass to the consideration of vibrations of orders 2, 3, ..., 
expecting (as, in fact, is found to be the case) that instability will first set in through 
a vibration of order 2. 

It is only possible to discuss special cases, and the one which is most amenable to 
analysis is that in which the pressure and density are connected by Laplace's law, 
p = c(p 2 —or 2 ). It is first proved (§§ 8-11) that, for a mass of such matter at rest, 
the spherical form is stable for all displacements. Later (§§ 15-22) it is shown that 
this is true when c varies inside the mass ; it is true even up to the case which is the 
most likely to be unstable, in which the matter in the interior is of negligible density 
and the main part of the mass is collected in a surface crust— a sort of astronomical 
soap-bubble. 

We proceed next to examine for what amount of rotation these figures will become 
unstable, treating first the case in which c is the same throughout the mass. 
Imagining c and <j to vary we can get a variety of types of mass. The surprising 
result is obtained (by something short of strict mathematical proof) that the figure 
which is the first to become unstable (as w 2 \2irp increases uniformly for them all) is 
the perfectly incompressible one — -gravitational instability appears to act in the 
unexpected direction, at any rate when the degree of rotation is measured by iv 2 /27rp, 
p being the mean density. As it was not possible to obtain strictly accurate figures 
in this case, the result was checked by considering the artificial, but physically 
analogous, problem of rotating cylinders of Laplacian matter, in which it was possible 
to obtain perfectly exact results (§ 14). The result was confirmed, and the additional 
information was obtained that the value of w 2 \2irp remains surprisingly steady 
through quite a wide range of compressibility (vide table on p. 471). 

The physical reason for this can, I think, be understood as follows. The more 
compressible the matter is the more it tends to concentrate near the centre, i.e., in 
just those regions where the " centrifugal force" obtains, so to speak, least grip on it. 
Incompressibility neutralises the gravitational tendency to instability, but tends to 
compel the matter to place itself so that the rotational tendency to instability can 
act at the best advantage. 

The similar problem is next investigated (§§ 23, 24) when c varies inside the mass ; 
in particular, the limiting case of a soap-bubble-like mass is considered. Again the 
surprising result emerges that the value of iv 2 /27rp needed to establish instability of 
the symmetrical configuration is just about the same as before (vide table, p. 481). 
The matter is now constrained to remain, so to speak, on the rim of a fly-wheel where 
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the centrifugal force can act at the best advantage and gravitational instability has 
full scope. If p is the mean density of the crust, w 2 /2irp must obviously be less than 
before. But if p is the mean density of the whole mass, p/p is also much smaller. 
These two quantities march with approximately equal steps, so that %v 2 J2irp remains 
almost unaltered. 

Thus we have the general result that for all the varied types of mass that have 
been considered the spheroidal or quasi-spheroidal form always becomes first unstable 
for just about the same value of w 2 J2irp. If, from the point of view of discovering 
iiew processes in nature, the present investigation has been somewhat barren, at least 
we may reflect that the work of Darwin and Poincake has been shown, to some 
extent, to have an enhanced value, in that it seems to apply to the seal bodies of 
nature and not merely to mathematical abstractions. 



